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Using the time-dependent density matrix renormalization group (TDMRG) we theoretically study 
the collision of a one-dimensional gas of interacting bosons with an impurity trapped in a shallow 
box potential. We study the dynamic response of both the impurity and the bosonic gas and find 
an approximate independence of the sign of the interaction (attractive or repulsive) between the 
impurity and the bosons. This sign-independence breaks down when the interaction energy is of the 
same order as the box potential depth, leading to resonant pair tunneling. Our predictions can be 
tested using contemporary techniques with ultracold atoms. 
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Introduction. — Ultracold atoms are an ever more ver¬ 
satile toolbox for simulating complex condensed matter 
systems as well as for experimentally realizing various 
theoretical “toy models.” In particular, there has been 
significant interest in one-dimensional (ID) systems over 
the previous decade. Such systems were previously re¬ 
garded mainly as a mathematical curiosity, but were 
realized experimentally in the early 2000s m- One¬ 
dimensional systems can be created by tightly confining 
a dilute cloud of atoms in two of the three dimensions 
0] and the strength of the inter-particle interactions can 
be tuned using Feshbach resonances [5] and confinement- 
induced resonances HE- Subsequent experimental ef¬ 
forts in ID include various realizations of Bose gases and 
the dynamics therein mm, studies of the dynamics of 
impurities [T7l - l20| . and fermionic few-body systems 
[23] . Theoretically, one-dimensional systems are of inter¬ 
est because of the existence of exact solutions and power¬ 
ful approximative techniques for various relevant many- 
body scenarios - see, e.g., Refs. [24lf27] for reviews. In 
addition, the ground state solutions of several classes of 
one-dimensional Hamiltonians and the dynamics of par¬ 
ticles described thereby can be calculated numerically es¬ 
sentially to arbitrary accuracy using the density matrix 
renormalization group [25], time-evolving block decima¬ 
tion [251 or exact diagonalization of the Hamiltonian m- 

The effect of an impurity on the properties of a ID 
quantum gas has been extensively studied [5TH41) ; in 
higher dimensions one can experimentally realize the so- 
called Fermi polaron 02] 03], described accurately using 
a simple variational model [44] . This quasiparticle de¬ 
scription works even in ID 05IH7] . even though strictly 
speaking these quasiparticles are not well-defined in ID. 
A bosonic analog, the Bose polaron , has also been pro¬ 
posed [48l ~ F5(7j and experimentally realized in ID (20]. 

Here, we will consider a scenario depicted schemati¬ 
cally in FIG. [T] An ensemble of interacting ID bosons 
is initially trapped in a harmonic potential (for related 
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FIG. 1. (color online). Schematic depiction of the problem 
under consideration. Initially, Nf majority component bosons 
(red dashed line and arrows) are trapped in a harmonic trap, 
while a single j, impurity (blue dashed line and arrow) is 
trapped in a separate box potential. At a time t = 0, the 
harmonic trap is instantaneously displaced (dotted line), al¬ 
lowing the jvparticles to collide and interact with the trapped 
impurity, x indicates the one-dimensional position. 


theoretical studies of ID bosons see Refs. [511152) ). At 
a time t = 0, the harmonic trap is displaced, allowing 
the bosons to collide and interact with an impurity. This 
entails a quantum quench (a sudden change in one of 
the parameters of the Hamiltonian); such quenches have 
been studied in particular in the context of bridging the 
gap between the micro-world of quantum physics and 
macroscopic statistical ensembles f53H55j . We place the 
impurity in a finite well or “box” potential m and in¬ 
vestigate the tunneling dynamics [23j 1681171) and spin 
transport [72] of the particle outside the box as a func¬ 
tion of the various parameters of the setup, such as the 
interaction between the majority bosons and the depth of 
the box potential. Recently, a similar system consisting 
of a few bosons interacting with an impurity on the sites 
of an optical lattice was studied experimentally [73]. 

Model. — We consider a one-dimensional system at zero 
temperature with an impurity (4) interacting with a num¬ 
ber of (t) bosons, which themselves can interact with 
each other. This system can be described using the Bose- 
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FIG. 2. (color online). Density profiles for the impurity nj. (top panels) and the majority component bosons n-|- (lower panels) 
and interaction strengths (columns) U/J = ±1,±10. In the panels, the s-axis indicates position and the y-axis time. Darker 
colors indicate higher densities. Note the many-body effects for U/J = ±1 as evidenced by the formation of “branches” in 
the density profile. In the case of a single majority component particle (right panels), only two such branches are visible. A 
“quantum Newton’s cradle”-type effect is visible in the majority component response [lT] (lower panels, first four columns). For 
strong interactions U/J = ±10, the impurity is fully kicked out of the box potential, unless it is sufficiently strongly trapped 
(fifth column). The number of majority component bosons Nf = 4, the intra-species interaction Uf/J = 5 and the depth of the 
box potential V\,/J = 0.05 unless otherwise noted. See also the online supplementary material for animated density profiles. 


Hubbard Hamiltonian for a two-component system: 

H = - J b\ i b cri+ 1 + H.c. + C/ t ^2 “ x ) 

cri i 

± bJ ^ ± ^ ^ Fc 77(72 (1) 

i rri 

Here J represents the hopping parameter, £/| is the on¬ 
site interaction energy between the majority component 
t-bosons, U represents the interaction between the im¬ 
purity and the f-bosons and fid = b' aj bd is the den¬ 
sity operator, where b^ is the bosonic annihilation (cre¬ 
ation) operator for a species of generalized spin er g {f 
, j,} at lattice site i. The external potential is added 
through an on-site term V al , which is different for both 
species. We consider only the case where the on-site 
density ( fi a )i <C 1 V«, in which case the Bose-Hubbard 
Hamiltonian can be mapped to its continuum equivalent, 
the (two-component) Lieb-Liniger gas with contact (£- 
function) interactions PH [71] with added terms for the 
external potentials. In this mapping, the hopping pa¬ 
rameter J is mapped to the mass m of a particle, which 
we assume to be the same for both species. The on-site 
interaction energies Uf and U can be associated with the 
strength of the intra- and inter-species contact interac¬ 
tions respectively. 

Results. — We use the time-dependent density ma¬ 
trix renormalization group (TDMRG) [28]. This method 
yields a numerically exact solution for the ground state 
and time dynamics of the Bose-Hubbard Hamiltonian 


([I]). Although we are interested primarily in the con¬ 
tinuum limit, we will state the Hubbard model pa¬ 
rameters for transparency. Our quench protocol is de¬ 
picted schematically in FIG. [T] We use a system of 
L = 250 sites, intra-species interactions U^/J = 5 and 
vary the inter-species interaction U and number of ma¬ 
jority component bosons N^. The choice CA/J = 5 
approximately corresponds to a Lieb-Liniger parameter 
7 ll = U^-/2Jh^ « 25, using half the peak density fq- 
of the majority bosons. Initially, the majority bosons 
are trapped in a harmonic trap V^/J = 10 _4 (sa — *) 2 
and the impurity is trapped in a box potential of depth 
Vj J J = 0.05 for 85 < i < 95. We prepare these sub¬ 
systems in their respective ground states with U = 0. 
To elucidate the effects of the interactions between the 
particles, we assume that only the impurity “feels” the 
trapping potential of the box potential, while only the 
majority bosons feel the harmonic trap. Such a setup 
might be realized experimentally using species-selective 
trapping potentials m- The harmonic trap is displaced 
from x\ = 50 to xb = 80 at a time t = 0, while at 
the same time the inter-species interaction U is switched 
on. This choice ensures that the initial overlap between 
the two clouds is small enough so that the effect of the 
interaction quench is not relevant. Furthermore, this 
particular quench protocol implies that, in the contin¬ 
uum limit and in the absence of interactions with the 
impurity, the cloud of majority component particles will 
perform harmonic motion about xb indefinitely. This 
allows us to see precisely the effect of the impurity on 
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the majority component particles. The impurity is ini¬ 
tially only weakly confined; the box potential permits 
exactly one single-particle bound state for V^/J = 0.05. 
After the displacement of the harmonic trap, we follow 
the dynamics of the system until t = 250 J -1 . This time 
corresponds approximately to the period of the oscilla¬ 
tion of the majority component bosons in the harmonic 
trap. The density profiles for the impurity n ^ = {h]/)i 
and the majority bosons rif = (h-f)i as a function of time 
for various parameters are depicted in FIG. [2j 

After the displacement of the harmonic trap, the cloud 
of the majority propagates towards the impurity, which 
then has a certain probability of being kicked out of the 
shallow box trap. We can quantify this probability by 
computing the following quantity: 

( 2 ) 

ij 

where is the reduced density matrix of the .(-particle 
and (j> is the single-particle ground state of the impurity 
obtained through an exact diagonalization method inde¬ 
pendent from our TDMRG implementation. This quan¬ 
tity is equal to 0 initially as the impurity is found in the 
ground state of the box potential, and reaches values of 
up to unity in the case the particle is completely kicked 
outside the box potential. In addition, we track the time 
evolution of the purity of the reduced density matrix: 

7 = Tr(p|). ( 3 ) 

This quantity enables us to quantify the degree of entan¬ 
glement between the majority bosons and the impurity. 
It is equal to unity at t = 0 since the system is prepared 
such that the impurity and the f-bosons are independent, 
i.e. the initial state is a pure state. The quantities © and 
^ are depicted in FIG. [3] 

Interestingly, the behavior of the probability V can be 
described remarkably well using a simple two-particle pic¬ 
ture. In free space, the reflection coefficient (the square 
of the scattering amplitude) \f(k)\ 2 for a particle with 
momentum k interacting through a contact potential 
u{\x — x'\) = gS(\x — x'\) with a second, stationary parti¬ 
cle at position x' is given by (we choose units 

where Ti = m/2 = 1): 

^ 2 = 7T¥- (4) 

We should therefore expect the tunneling of the impu¬ 
rity outside of the trap, i.e. the reflection of the impurity 
from the majority component particles, to scale as U 2 
for weak interactions (as might also be expected from 
a perturbative approach), with a saturation to complete 
reflection for \U/J\ —► oo. This behavior is recovered by 
the numerical results in the limits of weak and strong 
interaction. Interestingly, an interaction quench leads 
to similar results m- Note that this picture ignores the 
depth of the box potential confining the impurity; indeed, 


for a sufficiently deep potential (see FIG.[2j fifth column) 
tunneling of the impurity is suppressed, as expected. In 
this case there is also an additional single-particle bound 
state in the box potential. The effect of Vj, is examined 
more closely in FIG. ©>• Note that the effect of the sign 
of the interaction is more pronounced for Vj./ J = 0.2. In 
the weakly interacting limit, the U 2 scaling is recovered, 
albeit with a suppressed tunneling probability. Mean¬ 
while, in the strongly interacting limit U —» oo, the tun¬ 
neling probability saturates to a finite value smaller than 
unity. However, the sign independence remains in both 
the weakly and strongly interacting limits. 

For strongly attractive interactions, no significant pair¬ 
ing occurs between the f bosons and the impurity, and 
the j, particle is simply kicked outside the well with 
V approaching unity (corresponding to reflection), even 
though the static ground state of the system consists of 
the impurity pairing with one of the majority component 
bosons in the harmonic trap. This effect was observed in 
an experiment with spin-polarized fermionic clouds m- 
A similar picture was used to analyze the collision of 
fermion clouds in Ref. [77], although the sign of the in¬ 
teraction was not varied in this study. In Ref. [76], a 
semiclassical Boltzmann approach was used to analyze 
the collision of two spin-polarized fermionic clouds, also 
predicting an independence of the sign of the interaction, 
while the authors of Ref. m use a hydrodynamic ap¬ 
proach and draw similar conclusions. In lattice models, 
theoretical approaches mm and an experiment [82] 
yield similar effects. Note, however, that in the experi¬ 
ment of Ref. [82] and the associated theory, the interac¬ 
tion sign independence relies on a specific symmetry of 
the lattice dispersion relation that is not present in the 
situation we discuss here. 

In our simulations with spin-polarized bosons, we find 
the sign symmetry to hold over a wide range of interac¬ 
tions U/J. However, the independence of the sign does 
not hold universally. The most significant deviation oc¬ 
curs for intermediate interactions U/J « 1 [80]. This 
also manifests itself in the stronger entanglement of the 
impurity with the majority component bosons as shown 
in FIG. [sjc, as measured through the purity ©• We at¬ 
tribute this behavior to the relative importance of pair 
formation when the interaction energy is of the same 
order as the single-particle ground state energy E 0 in 
the box potential. To quantify this, consider the map¬ 
ping of the Bose-Hubbard Hamiltonian ([I]) to the Lieb- 
Liniger Hamiltonian. In this continuum mapping, we 
can associate the hopping J with the mass of a particle 
through J — > h 2 /2md (38], where d is the lattice spac¬ 
ing, while U —3 g/d, where g measures the strength of 
the ID contact interaction. The energy of a bound state 
dimer of two particles interacting through a delta func¬ 
tion potential is —mg 2 /Ah 2 = — U 2 /8J. Equating this 
energy to Eq, which we can obtain numerically, we ob¬ 
tain for Vj J J = 0.05 that E 0 /J ss —0.027, which yields 
U/J~ —0.46, in good agreement with FIG. §3. Mean¬ 
while, for Vj J J = 0.2 we obtain U/ J ~ —1.1, confirming 
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FIG. 3. (color online), a) probability V that the impurity 
is kicked outside the box potential, eq. @. b) (note the 
double log scale) V evaluated at a fixed time t = 200 J -1 for 
various interaction strengths U/J and box potential depths 
V\./J. Both positive (lower lines) and negative (upper lines) 
values of U/ J are shown. The dotted line is a guide to the 
eye and is proportional to I/ 2 . Note that for V\,/J = 0.1, 0.2 
we subtracted the overlap with the first two single-particle 
eigenstates, c) purity 7 of the reduced density matrix p 4 ., 
Eq. (§, as a function of time. The strongest entanglement 
of the impurity and the majority component bosons occurs 
for intermediate interactions U/J = ±0.3,1. Nf = 4, V 4 ./J = 
0.05 for all panels unless noted otherwise. 


the shift to higher interaction strengths for deeper box 
potentials. Thus, the independent collision model (ICM) 
of Ref. [77], employing a series of independent, consec¬ 
utive two-body collisions, breaks down when pairing ef¬ 
fects are important. For strong interactions \U/J\ 1 

the purity ([3]) saturates at a finite value close to unity. 
We explain this behavior through the impurity being 
kicked out of the box potential completely and then prop¬ 
agating as a free particle for later times. Hence, while the 
interaction between the majority bosons and the impu¬ 
rity is strong, the induced entanglement is not. 

The formed pairs are short-lived and dissociate due 
to the harmonic trapping that affects only the majority 
component bosons. Therefore, the bosons act as an in¬ 
termediary, enhancing the tunneling of the impurity out 
of the trap as is visible in the first two panels of the top 
row in FIG. [2] We can visualize this effect by considering 


the pair density ffl] : 

V = V n t n -1 “ n !T> ( 5 ) 

where the square root of the doublon density = 
yj(b^pib^b-f), and n a = {b\b a ) as before. When 77 < 0, 
the density of doublons is greater than the product of 
densities, which indicates pairing. Conversely, when 
> 0 we have anti-bunching. We consider the case where 
we expect the largest difference between repulsive and at¬ 
tractive interactions based on FIG. §>> i.e. U/J = ±0.3. 
In FIG. [4] (lower panel) the process of pairing followed by 
dissociation is clearly visible. Initially, when the bosonic 
cloud collides with the impurity, pairs are formed as in¬ 
dicated by the negative rj. The pair propagates briefly 
before dissociating, which coincides with positive 77 at 
later times. 



FIG. 4. (color online). Density of the impurity (top panel), 
the majority component bosons n-f (middle panel) and the 
pair density r/ (lower panel) as a function of time and po¬ 
sition. Negative values in the lower panel indicate pairing, 
while positive values indicate anti-bunching. U/J = —0.3, 
Uf/J = 5, Vj J J = 0.05 and Nf = 4. 

The reflected branch of the impurity (see FIG. [2]) 
propagates as an undisturbed, expanding Gaussian wave 
packet with a constant propagation velocity of the center 
of mass. Intuitively, one would expect this propagation 
velocity to depend on the “kick” given to the majority 
bosons (xb — %a)- This quantity can be associated with 
the momentum k in Eq. Q. Our numerics appear to sup¬ 
port this picture for small variations around our choice 
of (o;b — £a) 5 although we cannot perform a systematic 
analysis of the dependence on {x-q — xa) due to numerical 
constraints. In addition, we find that for a number of ma¬ 
jority bosons Nf = 4 compared to a single f-boson, the 
impurity receives a larger kick, indicated by the steeper 
slope in FIG. [2j as it interacts with multiple particles. 
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FIG. 5. (color online). As in FIG. [4] but for U/J = 0.3. 


Conclusions .— We study the collision of a one¬ 
dimensional Bose gas with a weakly trapped impurity. 
The probability of the impurity being kicked out of the 
box potential depends primarily on the magnitude of the 
interaction between the majority component bosons and 
the impurity, but only weakly on the sign of the interac¬ 
tion. This is explained in terms of a two-particle scatter¬ 
ing picture. 

We predict a pair tunneling effect when the interac¬ 
tion energy and the single-particle ground state energy 
Eq of the impurity are of the same order. This can be 
compared to pair tunneling of pairs of attractively in¬ 
teracting distinguishable particles that are prepared in a 
trap [22 ESI IS] • In that case, the interaction leads to a 
reduced tunneling probability of the pair from the trap, 
as the (negative) interaction energy leads to an effective 


higher barrier height. Conversely, here we find an in¬ 
creased tunneling probability mediated by the attractive 
interaction because the energy released when the pair 
forms is converted to kinetic energy of the pair, allow¬ 
ing it to tunnel outside the box potential. This feature 
is broad (see FIG. [3 }d) because the majority component 
bosons have a large spread in kinetic energy. Further¬ 
more, our results indicate that the enhanced tunneling 
peak moves towards higher interactions for deeper box 
potentials (lower Eq ), consistent with the pair tunneling 
picture. 

The majority component bosons are found near the 
Tonks-Girardeau limit U^/J —> oo. In this limit, 
“fermionization” occurs, where local observables are 
identical to a non-interacting spinless gas of fermions. 
However, the momentum distribution of the Tonks- 
Girardeau gas is very different, exhibiting a characteristic 
1 /g 1 2 3 4 5 6 7 8 decay, where q is momentum 1811 . It would be of in¬ 
terest to see whether this affects the problem considered 
here. 

One aspect that is missing in the picture of Eq. Q is 
the depth of the box potential V±. Indeed, for the impu¬ 
rity to tunnel out of the box potential, it needs an addi¬ 
tional energy ss Eq. This opens up the prospect of using 
the scheme outlined here to measure the off-shell scatter¬ 
ing amplitude (proportional to the two-body f-matrix), 
where the energy is not conserved in the collision. 
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